On the structures in the afterglow peak emission of gamma ray 

bursts 



Remo Ruffini^'^, Carlo Luciano Bianco^'^, Pascal Cliardonnet^'^, Federico Fraschetti^'^, 

She-Slieng Xue^'^ 

Vn" 

o : abstract 
o 

(N 



> 

oo 



Using GRB 991216 as a prototype, it is shown that the intensity substructures 
Q ■ observed in what is generally called the "prompt emission" in gamma ray bursts 

O ' (GRBs) do originate in the collision between the accelerated baryonic matter 

^ ■ (ABM) pulse with inhomogeneities in the interstellar medium (ISM). The initial 

phase of such process occurs at a Lorentz factor 7 ~ 310. The crossing of ISM 
inhomogeneities of sizes AR ~ 10^^ cm occurs in a detector arrival time interval 
^ ■ of ~ 0.4 s implying an apparent superluminal behavior of ~ lO^c. The long 

Q ■ lasting debate between the validity of the external shock model vs. the internal 

^ ' shock model for GRBs is solved in favor of the first. 

O '■ 

^ ■ Subject headings: black hole physics — gamma rays: bursts — gamma rays: 

I ' observations — gamma rays: theory — ISM: clouds — ISM: structure 

|: 

! To reproduce the observed light curve of GRB 991216, we have adopted, as initial 

j> ! conditions (Ruffini et al. 2002a) at t = 10~^^ s ~ s, a spherical shell of electron-positron- 
^ ! photon neutral plasma laying between the radii tq = 6.03 x 10^ cm and ri = 2.35 x 10^ cm: 
^ ! the temperature of such a plasma is 2.2 MeV, the total energy Etot = 4.83 x 10^^ erg and the 
total number of pairs N^+e- = 1.99 x 10^*^. 

Such initial conditions follow from the EMBH theory we have recently developed based 
on energy extraction from a black hole endowed with electromagnetic structure (EMBH) 
(Ruffini 1998; Preparata et al. 1998; Ruffini et al. 1999, 2000; Bianco et al. 2001; Ruffini et al. 
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2001a,b,c, 2002a), ro being the horizon radius, ri the dyadosphere radius and E'tot coinciding 
with the dyadosphere energy E^ya- The above set of parameters is uniquely determined by 
the value of E^ya- The EMBH energy (Christodoulou & Ruffini 1971) is carried away by 
a plasma of electron-positron pairs created by the vacuum polarization process (Damour & 
Ruffini 1975) occurring during the gravitational collapse leading to the EMBH (Cherubini 
et al. 2002; Ruffini & Vitagliano 2002). Such an optically thick electron-positron plasma self 
propels itself outward reaching ultrarelativistic velocities (Ruffini et al. 1999), interacts with 
the remnant of the progenitor star and by further expansion becomes optically thin (Ruffini 
et al. 2000). The physical reason for such an extraordinary process of self- acceleration, 
achieving in a tenth of seconds in arrival time an increase in the Lorentz gamma factor from 
7 = 1 to 7 ~ 300. has been shown to be critically dependent on E^ya and on the amount 
of baryonic matter engulfed by the plasma in its expansion (see Ruffini et al. 1999, 2000). 
It is interesting that this process is extremely efficient even in the present case, regardless 
of the relatively slow random thermal motion of the 2.2 MeV e+e^ plasma (see Ruffini et 
al. 2002a). As the transparency condition is reached, a proper GRB (P-GRB) is emitted 
as well as an extremely relativistic shell of accelerated baryonic matter (ABM pulse). It is 
this ABM pulse which, interacting with the interstellar medium (ISM), gives origin to the 
afterglow (see Ruffini et al. 2001b, 2002a). 

One of the most novel results of the EMBH model has been the identification of what 
is generally called the "prompt emission" (see e.g. Piran 1999, and references therein) as 
an integral part of the afterglow: the extended afterglow peak emission (E-APE) (Ruffini 
et al. 2001a,b, 2002a). This result is clearly at variance with the models explaining the 
"prompt emission" with ad-hoc mechanisms distinct from the afterglow process (see e.g. 
Rees & Meszaros 1994; Meszaros & Rees 1997; Kobayashi et al. 1997; Rees & Meszaros 1998; 
Meszaros & Rees 2001; Rees & Meszaros 2000; Kumar & Piran 2000; Meszaros 2002). The 
fact that the EMBH model, using GRB 991216 as a prototype, has allowed the computation 
of the temporal separation of the P-GRB and the E-APE to an accuracy of a few milliseconds 
and also to predict their relative intensities within a few percent can certainly be considered 
a major success of the model (see Ruffini et al. 2001a,b, 2002a). 

The aim of this letter is to report a further extension of the EMBH model in order 
to identify the physical processes giving origin to the intensity variability observed in the 
El- APE on time scales as short as a fraction of a second (Fishman & Meegan 1995), which 
contrasts with the smoother emission in the last phases of the afterglow (see e.g. Costa et 
al. 2001). 

In our former work on the EMBH model (Ruffini et al. 2001b, 2002a), we have assumed 
an homogeneous ISM with a density riism —< nism >— 1 particle/cm^ and we have also 
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assumed that during the coUision of the ABM pulse with the ISM the "fully radiative con- 
dition" applies. These assumptions have led to the theoretical prediction of the power-law 
index of the afterglow slope n = —1.6 in excellent agreement with the observational data 
n = —1.616 ± 0.067 (Halpcrn ct al. 2000). Our goal here is to show that the variability in 
the E-APE can indeed be traced back to inhomogeneities in the ISM. We again consider, 
like in the previous work, the case of an ABM pulse expanding with spherical symmetry 
(i.e. no beaming) and for simplicity we describe the ISM inhomogeneities as spherical shells 
concentric to the ABM pulse. Each shell has a selected density and a constant thickness 
= 1.0 X 10^^ cm. 

We recall now the relation between the relativistic beaming angle and the arrival time 
of the emitted photon on the detector. The visible part of the ABM pulse spherical surface 
is constrained by: 

COS'&>-, (1) 
c 

where is the angle in the laboratory frame between the radial direction of each point on 
the ABM pulse surface and the line of sight and v is the expansion speed (Ruffini et al. 
2002b). This follows from the requirement that in the comoving frame the component of 
the photon momentum along the radial expansion velocity direction be positive, in order to 
escape. There exists then a maximum allowed d value "dmax defined by cos'i^^a^; = {v/c) (see 
Fig. la). 

Due to the high value of the Lorentz 7 factor (~ 310) for the bulk motion of the ABM 
pulse, the spherical waves emitted from its external surface do appear extremely distorted to 
a distant observer. To show this we need to express the photon arrival time at the detector 
of a function of its emission time t and angle i). We set t = when the plasma starts to 
expand, so that r (0) = r^s- We then have (see Ruffini et al. 2002b): 

tt = {l + z) ft-£^ii^>^cos^+^ 

" y c c 

where z is the redshift of the source. Then, in order to compute the arrival time of the 
emitted radiation, we must know all the previous values of the source velocity starting from 
t = 0. The great advantage of the EMBH model is that for the first time we have been 
able to obtain the precise values of the gamma Lorentz factor as a function of the radial 
coordinate or cquivalcntly of the laboratory time (sec Fig. 2). This allows us, for the first 
time, to evaluate Eq.(2) and correspondingly determine the surfaces that emits the photons 
detected at a fixed arrival time f^, which we will call in the following "equitemporal surfaces" 
(EQTS). The profiles of such surfaces are reported in Fig. lb. We emphasize once again the 
direct connection between the evaluation of the EQTS and the entire past history of the 
source. 
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We have created an ISM inhomogeneity "mask" (see Fig. 3 and Tab. 1) with the main 
criteria that the density inhomogeneities and their spatial distribution still fuUfiU < riism >— 
1 particle/cm^. 

The source luminosity in a detector arrival time and per unit solid angle dO, is given 
by (details in Ruffini et al. 2002b): 



where Ae is the energy density released in the interaction of the ABM pulse with the ISM 
inhomogeneities measured in the comoving frame, A = 7(1 — {v/cjcosi}) is the Doppler 
factor and rfS is the surface element of the EQTS at detector arrival time on which the 
integration is performed. In the present case the Doppler factor A^'^ in Eq. (3) enhances the 
apparent luminosity of the burst, as compared to the intrinsic luminosity, by a factor which 
at the E-APE is in the range between 10^° and 10^^! 

The results are given in Fig. 5. We obtain, in perfect agreement with the observations 
(see Fig. 4): 

1. the theoretically computed intensity of the A, B, C peaks as a function of the ISM 
inhomogneities; 

2. the fast rise and exponential decay shape for each peak; 

3. a continuous and smooth emission between the peaks. 

Interestingly, the signals from shells E and F, which have a density inhomogeneity 
comparable to A, are undetectable. The reason is due to a variety of relativistic effects and 
partly to the spreading in the arrival time, which for A, corresponding to 7 = 303.8 is 0.4s 
while for E (F) corresponding to 7 = 57.23 (56.24) is of 10.2 s (10.6 s) (see Tab. 1 and Ruffini 
et al. (2002b)). 

In the case of D, the agreement with the arrival time is reached, but we do not obtain the 
double peaked structure. The ABM pulse visible area diameter at the moment of interaction 
with the D shell is ~ 1.0 x 10^^ cm, equal to the extension of the ISM shell (see Tab. 1 and 
Ruffini et al. 2002b). Under these conditions, the concentric shell approximation does not 
hold anymore: the disagreement with the observations simply makes manifest the need for 
a more detailed description of the three dimensional nature of the ISM cloud. 

The physical reasons for these results can be simply summarized: we can distinguish two 
different regimes corresponding in the afterglow of GRB 991216 respectively to 7 > 150 and 
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to 7 < 150. For different sources tfiis value may be sligfitly different. In tfie E-APE region 
(7 > 150) tlie GRB substructure intensities indeed correlate with the ISM inhomogcneities. 
In this limited region (see peaks A, B, C) the Lorentz gamma factor of the ABM pulse 
ranges from 7 ~ 304 to 7 ~ 200. The boundary of the visible region is smaller than the 
thickness AR of the inhomogeneities (see Fig. 1 and Tab. 1). Under this condition the 
adopted spherical approximation is not only mathematically simpler but also fully justified. 
The angular spreading is not strong enough to wipe out the signal from the inhomogeneity 
spike. 

As wc descend in the afterglow (7 < 150), the Lorentz gamma factor decreases markedly 
and in the border line case of peak D 7 ~ 140. For the peaks E and F we have 7 ~ 50 and, 
under these circumstances, the boundary of the visible region becomes much larger than the 
thickness AR of the inhomogeneities (see Fig. 1 and Tab. 1). A three dimensional description 
would be necessary, breaking the spherical symmetry and making the computation more 
difficult. However we do not need to perform this more complex analysis for peaks E and F: 
any three dimensional description would a fortiori augment the smoothing of the observed 
flux. The spherically symmetric description of the inhomogeneities is already enough to 
prove the overwhelming effect of the angular spreading (Ruffini et al. 2002b). 

On this general issue of the possible explanation of the observed substructures with the 
ISM inhomogeneities, there exists in the literature two extreme points of view: the one by 
Fenimore and collaborators (see e.g. Fenimore et al. 1996, 1999; Fenimore 1999) and Piran 
and collaborators (see e.g. Sari & Piran 1997; Piran 1999, 2000, 2001) on one side and the 
one by Dermer and collaborators (Dermer 1998; Dermer et al. 1999; Dermer & Mitman 1999) 
on the other. 

Fenimore and collaborators have emphasized the relevance of a speciflc signature to be 
expected in the collision of a relativistic expanding shell with the ISM, what they call a fast 
rise and exponential decay (FRED) shape. This feature is confirmed by our analysis (sec 
peaks A, B, C in Fig. 5). However they also conclude, sharing the opinion by Piran and 
collaborators, that the variability observed in GRBs is inconsistent with causally connected 
variations in a single, symmetric, relativistic shell interacting with the ambient material 
("external shocks") (Fenimore et al. 1999). In their opinion the solution of the short time 
variability has to be envisioned within the protracted activity of an unspecified "inner engine" 
(Sari & Piran 1997); see as well Rees & Meszaros (1994); Panaitescu & Meszaros (1998); 
Meszaros & Rees (2000, 2001); Meszaros (2002). 

On the other hand, Dermer and collaborators, by considering an idealized process occur- 
ring at a fixed 7 = 300, have reached the opposite conclusions and they purport that GRB 
light curves are tomographic images of the density distributions of the medium surrounding 
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the sources of GRBs (Dermer & Mitman 1999). 

Prom our analysis we can conclude that Dermer's conclusions are correct for 7 ~ 300 
and do indeed hold for 7 > 150. However, as the gamma factor drops from 7 ~ 150 to 
7 ~ 1 (see Fig 2), the intensity due to the inhomogeneities markedly decreases also due to 
the angular spreading (events E and F). The initial Lorentz factor of the ABM pulse 7 ~ 310 

decreases very rapidly to 7 ~ 150 as soon as a fraction of a typical ISM cloud is engulfed (see 
Fig. 2 and Tab. 1). We conclude that the "tomography" is indeed effective, but uniquely in 
the first ISM region close to the source and for GRBs with 7 > 150. 

One of the most striking feature in our analysis is clearly represented by the fact that 
the inhomogeneities of a mask of radial dimension of the order of 10^^ cm give rise to arrival 
time signals of the order of 20 s. This outstanding result implies an apparent "superluminal 
velocity" of ~ lO^c (see Tab. 1). The "superluminal velocity" here considered, first intro- 
duced in Ruffini et al. (2001a), refers to the motion along the line of sight. This effect is 
proportional to 7^. It is much larger than the one usually considered in the literature, within 
the context of radio sources and microquasars (see e.g. Mirabel & Rodriguez 1994), referring 
to the component of the velocity at right angles to the fine of sight (see details in Ruffini et 
al. 2002b). This second effect is in fact proportional to 7 (see Rees 1966). We recall that this 
"superluminal velocty" was the starting point for the enunciation of the RSTT paradigm 
(Ruffini et al. 2001a), emphasizing the need of the knowledge of the entire past worldlines of 
the source. This need has been ftirther clarified here in the determination of the EQTS sur- 
faces (see Fig. lb) which indeed depend on an integral of the Lorentz gamma factor extended 
over the entire past worldlines of the source. In turn, therefore, the agreement between the 
observed structures and the theoretical predicted ones (see Figs. 4-5) is also an extremely 
stringent additional test on the values of the Lorentz gamma factor determined as a function 
of the radial coordinate within the EMBH theory (see Fig. 2). 

We thank M. Rees for pointing out the necessity of presenting our results, R. Giacconi 
for suggestions on the wording of the manuscript and an anonymous referee for excellent 
advices. 
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Fig. 1. — a) We represent the temporal evolution of the ABM pulse visible area. The dashed 
half-circles are the expanding ABM pulse at selected values of its radius corresponding to 
different laboratory times. The black curve marks the boundary of the visible region. The 
EMBH is located at position (0,0). The X (Y) axis is directed along (orthogonal to) the hne 
of sight. In the earliest GRB phases the visible region is squeezed on the hne of sight, while 
in the latest afterglow phases almost all the emitted photons reaches the observer, b) In 
the same coordinate system used in a), we represent (solid lines) the equitemporal surfaces 
(EQTS) (sec text). They corresponds to values of the arrival time ranging from 5 s (the 
smallest surface on the left of the plot) to 60 s (the largest one on the right) in steps of 5 s. 
The dashed lines are the boundaries of the ABM pulse visible area. Note the different scale 
on the two axis, indicating the very high EQTS "effective eccentricity". The arrival time 
interval has been chosen to encompass the E-APE emission, occurring between ~ 15 s and 
~ 40 s (see Figs. 4-5 and Tab. 1). 
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Distance from the EMBH (cm) 



Fig. 2. — The theoretically computed gamma factors of the expanding pulse are given as a 
function of its distance from the EMBH. The dashed line corresponds to riism = 1 while the 
solid line corresponds to the density profile given in Fig. 3. The roman numerals corresponds 
to the different eras of the EMBH theory (see Ruffini et al. 2002a) . Near the E-APE (namely 
around Pa and Pl) the two curves differs markedly, due to the impact on the high density 
ISM regions which brakes the ABM pulse more efficiently. When the ABM pulse overcomes 
the ISM cloud the two curves coincide again, since the average density of the cloud is ~ 
1 particle/cm^. 
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Fig. 3.— The density profile ("mask") of an ISM cloud used to reproduce the GRB 991216 
temporal structure. As before, the radial coordinate is measured from the black hole. In 
this cloud we have six "spikes" with overdensity separated by low density regions. Each 
spike has the same spatial extension of 10^^ cm. The cloud average density is < riism >— 
1 particle/cm^. 
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Fig. 4.— The BATSE data on the E-APE of GRB 991216 (source: BATSE GRB hght 
curves 1999) together with an enlargement of the P-GRB data (source: BATSE Rapid Burst 
Response 1999). For convenience each E-APE peak has been labeled by a different uppercase 
Latin letter. 
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Fig. 5. — The source luminosity connected to the mask in Fig. 3 is given as a function of the 
detector arrival time (solid "spiky" line) with the corresponding curve for the case of constant 
nism — 1 particle/cm^ (dashed smooth hne) and the BATSE noise level (dotted horizontal 
hne). The "noise" observed in the theoretical curves is due to the discretization process 
adopted, described in Ruffini et al. (2002b), for the description of the angular spreading of 
the scattered radiation. For each fixed value of the laboratory time we have summed 500 
different contributions from different angles. The integration of the equation of motion of 
this system is performed in 22, 314, 500 contributions to be considered. An increase in the 
number of steps and in the precision of the numerical computation would lead to a smoother 
curve. 
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Table 1. For each ISM density peak represented in Fig. 3 we give the initial radius r, the 
corresponding comoving time r, laboratory time t, arrival time at the detector , diameter 
of the ABM pulse visible area d,,. Lorentz factor 7 and observed duration At^ of the 
afterglow luminosity peaks generated by each density peak. In the last column, the 
apparent motion in the radial coordinate, evaluated in the arrival time at the detector, 
leads to an enormous "superluminal" behavior, up to 9.5 x 10^ c. 



Peak 


r (cm) 


r(s) 


t{s) 




dv (cm) 


Atf (s) 


7 


"Superluminal" 


A 


4.50 X 10^^ 


4.88 X 10^ 


1.50 X 106 


15.8 


2.95 X 10^'' 


0.400 


303.8 


9.5 X lO-^c 


B 


5.20 X lO^^^ 


5.74 X 10^ 


1.73 X 10^ 


19.0 


3.89 X 10^4 


0.622 


265.4 


9.1 X lO^'c 


C 


5.70 X 10^^ 


6.54 X 10^ 


1.90 X 10^ 


22.9 


5.83 X 10" 


1.13 


200.5 


8.3 X 10*c 


D 


6.20 X 10^^ 


7.64 X 10^ 


2.07 X 10^ 


30.1 


9.03 X 10" 


5.16 


139.9 


6.9 X 10*c 


E 


6.50 X lO^*^ 


9.22 X 10^ 


2.17 X 10*5 


55.9 


2.27 X 10^5 


10.2 


57.23 


3.9 X lO^c 


F 


6.80 X 10^^ 


1.10 X 10'' 


2.27 X 10*3 


87.4 


2.42 X 10^5 


10.6 


56.24 


2.6 X lO^'c 



